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Abstract 

We propose a first order formalism for multi-centered black holes with flat tree- 
dimensional base-space, within the stu model of N = 2, D = 4 ungauged Maxwell-Einstein 
supergravity. This provides a unified description of first order flows of this universal sec- 
tor of all models with a symmetric scalar manifold which can be obtained by dimensional 
reduction from five dimensions. 

We develop a D = 3 Cartesian formalism which suitably extends the definition of 
central and matter charges, as well as of black hole effective potential and first order 
"fake" superpotential, in order to deal with not necessarily axisimmetric solutions, and 
thus with multi-centered and/or (under-)rotating extremal black holes. 

We derive general first order flow equations for composite non-BPS and almost BPS 
classes, and we analyze some of their solutions, retrieving various single-centered (static 
or under-rotating) and multi-centered known systems. 

As in the i 3 model, the almost BPS class turns out to split into two general branches, 
and the well known almost BPS system is shown to be a particular solution of the second 
branch. 



1 Introduction 



Multi-centered extremal black hole (BH) solutions in four-dimensional supergravity theories 
have been widely investigated in recent years [U El |3l HJ El El [7] . 

Motivated by the issue of matching the true BPS spectrum and the spectrum of spherically 
symmetric BHs in supergravity, in [8] a supersymmetric class of multi-centered BHs was in- 
troduced, whose BPS first order flow equations were then solved and analyzed in [9], revealing 
interesting features, such as fixed distance among the centers (in presence of mutually non-local 
electric- magnetic charge vectors). 

However, only in the last two years a group-theoretical approach [JJ, [2], based on nilpo- 
tent orbits and timelike reduction to three dimensions, allowed for a systematic construction 
and investigation of whole new classes of solutions, in which some or all BH centers are non- 
supersymmetric (non-BPS). 

An elegant approach to the flow dynamics of scalar fields in the background of single-centered 
extremal BH solutions of Maxwell-Einstein theories of (super) gravity, essentially based on the 
first order reformulation of the scalar equations of motion, was introduced in [10], and then 
developed in various works [TT1 IT2T [13]. The possibility to switch from second order to first 
order differential equations of motion - without doubling their number - has an applicative 
relevance. Indeed, due to the interplay between auxiliary fields and scalar charges, the first 
order formalism automatically discards blowing-up solutions; furthermore, the integration of 
first order equations is surely more manageable, and explicit forms of attractor flows can be 
more easily determined. 

The extension of such a formalism to non-supersymmetric multi-centered configurations was 
started in [7J by reducing the relevant action to a sum of squares, however without yielding 
explicit expressions for the flow equations and their corresponding governing functions. A 
consistent and explicit determination of the most general first order flow equations for non- 
BPS multi-centered and/or rotating BHs with flat three-dimensional base-space was achieved 
in [TJ], within the simplest model of iV = 2, D = 4 Maxwell-Einstein ungauged supergravity 
with a cubic prepotential, namely the so-called t 3 model, exhibiting only one vector multiplet, 
and whose uplift to D = 5 is "pure" minimal supergravity (see e.g. [T5]). 

In the present investigation, we further develop the approach of [JJ], and determine the 
general first order flow equations for all classes of multi-centered BH solutions with flat three- 
dimensional base-space in the so-called stu model [16]; in this model, three vector multiplets 
are coupled, in a triality-invariant way, to the N = 2 gravity multiplet, and the resulting com- 
pletely factorized rank-3 symmetric special Kahler manifold [SU(1, 1)/U(1)] 3 can be considered 
a universal sector of all symmetric scalar manifolds of D = 4, N ^ 2-extended supergravity 
theories which admit aD = 5 uplift. 

A key result is the reformulation of second order equations of motion in a manifestly D = 3 
Cartesian formalism, based on a timelike Lagrangian reduction D = 4 — > 3 in a stationary BH 
background [17], later specialized for flat spatial slices. This formalism allows for a consistent 
generalization of the BH effective potential [18] , and of its expression in terms of supersymmetry 
central charges and matter charges, as well as in terms of a first order "fake" superpotential 
[TO] , in not necessarily axisymmetric contexts, which thus can include multi-centered solutions. 

Various under-rotating stationary (BPS and non-BPS) single-centered solutions, as well as 
the known classes of BPS [8], composite non-BPS [JJ and almost BPS [JJ [3] multi-centered 
solutions (possibly constrained in a suitable way) are retrieved as particular solutions of the 
various first order systems which we propose. It should be stressed that, analogously to [14J, 
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in this paper we focus on the general formulation of first order systems of equations, leaving 
the determination of new sets of solutions and the investigation of their physical properties to 
future investigations. 

The paper is organized as follows. 

In Sec. 12], we start with the timelike reduction of Maxwell-Einstein-scalar action in the 
background of stationary metrics (then assumed to have flat spatial slices). As mentioned, this 
yields to a D = 3 Cartesian formalism, which is especially suitable to handle multi-centered 
and under-rotating single-centered solutions. 

Our general approach to first order formalism is then exploited in Sec. [3] for the various 
classes of multi-centered solutions, namely the BPS (Sec. I3.1.ip . almost BPS (Sec. I3.1.2p and 
composite non-BPS (Sec. 13. 1 .3[) classes, also retrieving and discussing known multi- and single- 
centered solutions. In particular, as found in [T4] in the t 3 model, we find that the almost BPS 
class splits into two branches, which can essentially be related to the BPS or non-BPS nature 
of the corresponding single-centered limits, as discussed in Sees. 15.11 and 15.21 

Finally, Sec. E] contains a summary and an outlook of results. 

Some notation and useful formulae are given in App. |A] 



2 Second Order Equations of Motion 

Let us consider the D = 4 Einstein-Maxwell action 



S = J d 4 x\/ — g 



~R + G a - a g^d,z a d u T + 1 vaxF^^ + I 



(2.1) 



which, depending on properties of the target space metric G a ^(z, z) and of the coupling ma- 
trices ijlae(z,z) and vaz,(z,z), may be the bosonic sector of some A" ^ 2-extended, D = 4 
Maxwell-Einstein supergravity. In the present investigation, we are interested in stationary so- 
lutions to the Einstein-Maxwell equations of motion; we will thus perform a timelike reduction 
to tree dimensions d la Breitenlohner- Gibbons- Mais on [17] . giving rise to a D = 3 Cartesian 
formalism, which is not necessarily axisymmetric. 

We start from a space-time metric satisfying the following Ansatz = 1, 2, 3): 

ds 2 = g^dx»dx u = e 2U(x \dt + Ui{x)dx 1 ) 2 - e~ 2U ^'y ij (x)dx i dx i (2.2) 

In this background, the Abelian vector potential A^ (defining the two-form field strength^ 
F^ u = 2d[f Ji A^) splits into D = 3 vector potentials af and Kaluza-Klein scalars b A : 

A^dx^ = b A {dt + Uidx 1 ) + afdx\ 
By using the D = 3 vectors af and Ui, one can construct the corresponding field strengths 

fij = 2<%4 Wij = 2d {i u, h 
respectively expressing the magnetic field and the rotation in tree dimensions. 



1 Throughout the paper, the (anti)symmetrization is defined with a 1/n! normalization if applied to n indices. 
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On stationary background (12.21) . the action (12. ip is equivalent to the following D = 3 one: 



S 



\ R (3) + G a - a ^d lZ a d 3 z a + Y j d t Ud 3 U - I e iU W l3 W^ 



+ 5 e-^^d^djb* - t e ^ AS (/A + h K Wij) ( f w + h * W i^ (2 .3) 
-VAxd t b A (7 Si + 6 s *W i )], 

where the D = 3 Hodge operator "*" is defined with respect to the 3-dimensional metric 7^: 

1 



"A 1 



z ijk A 



V7 



The D = 3 two-forms and Wjj can be dualized into scalars 6a and ip by adding the following 
terms to the action (12. 3p 



d x 



thus allowing for the two-forms to be expressed in terms of their duals: 

/Ay + b A w ij = ; l _ e -2U ^ £ m ^ fcfes + u ^,d k b^ ; 



W %3 



- iU e» k {d k ^ - b A d k b A ) 



(b) 



(2.4) 



as well as the action (12.3j) to be recast as: 
S= ( d 3 xy^ 



+ e~ 4 V {dtf - 6 A ^6 A ) (fi^ - b s djb s )] . 



(2.5) 



By grouping the Kaluza-Klein scalars b A and the scalars b A (dual to M) into a symplectic 
vector b a (with a running over contravariant and covariant symplectic index A) and introducing 
a scalar ijj 

b a =f b h \ ^ = 2^- b A b A 
the action (12.51) can be rewritten in more compact form: 

S I d 3 x\[^ 



(2.6) 



(2.7) 



For the definitions of the symplectic symmetric matrix M a p and the skew-symmetric sym- 
plectic product (■,■), see App. [A] Then, by introducing a symplectic vector /S = //vy) 



constructed from the field strengths, Eqs. ( 12 .4p enjoy a manifest symplectic covarianceH: 



fij = V7 e ilk ^'e-^M^n^dytf - b a W ij; (a) 

(b) 



(2.1 



The second order equations of motion pertaining to the action ( 12.71) consist of the Einstein 
equations 



= G a -ad {i z a d 3) -z a + diUdjU + I e~ 2U d l b a M al3 d j l 



+ r 



-41/ 



(2.9) 



and of the following ones: 



0, 

0; 



0. 



7 yie- w n a pbP djip + (b, djb) 



0; 



Tf %3 djU 



1 



--e-^^di^M^djb? 
- \ e" 4 u j ij (d^ + (b, d.b)) (d$ + (b, djb)) , 



(2.10) 



dG 



dM, 



Within the class of the stationary metrics ( 12. 2p . we further select those with D = 3 Euclidean 
metric 7^ = 5if 

(2.11) 



ds = e (dt + Uidx 



i\2 



-2U 



dx l dx\ 



Thus, we will henceforth restrict to consider D = 4 under-rotating BHs, with flat D = 3 
base space. When replacing (12. 2p with (12. lip . Einstein equations ( 12. 9 p become a first order 
constraint: 



Ga-ad^O^T + diUdjU + X - e-^diVMdjb + i e~ AU (d^ + (6, («9^ + (6, fyft) 



= 0, 

(2.12) 



and the equations of motion (12.101) get slightly simplified: 



di 



-4U 



d$ + (b, dib) 



0. 



$ [e~ 2U M^d^] =di e~ m n a pbP (d^ + (b,dib) 
dM - 



-\ e-^d^MdJ - 1 e- 4f/ (9^ + (b, d t b) 
di [Ga-ad^] = ^ d lZ %z l + \ e~ 2U dib^ d,b. 



(a) 
(b) 
(c) 
(d) 



(2.13) 



As a notation, let us introduce a D = 3 vector which, by virtue of (12.13a ). is divergenceless: 

Xl = e- 4C/ (d^ + (6, 9,6)) , <9 iX i = 0. (2.14) 



2 Eq. (|23Ji) is just the upper half of Eq. ([2 



Maxwell equations (12. 13b ) can immediately be integrated, yielding: 

dib a = e 2U M aP n prY H7, H« = H« + b a Xi, 0^ = 0, (2.15) 

where another D = 3 divergenceless vector Hf has been defined. The electromagnetic two-form 
field strength ff- and the two- form Wij corresponding to rotation are then just Hodge duals of 

the vectors Xi and #f defined in ([2TH)) - (j2TT5] ): 

= —EijkXk, fij = —£ijkH£ — b a Wij = —EijkH^- 

This allows one to rewrite the Einstein equations (12.121) as well as the rest of the equations of 
motion (12. 13ft as follows: 



G a - a d {l z a d 3) z & + diUdjU - %e 2U + \ e w X iXj = 0, (2.16) 

(2.17) 



2 ^ Aj) 



di (Gaad^) = d a G hl d lZ b da b + e 2U d a V n . 
In fl2TTB>(l2T7]) the D = 3 Cartesian tensor 

% = —HJMHj (2.18) 

has been introduced; it will be referred to as black hole potential. In fact, it generalizes, for the 
whole class of metrics (12.111) . the well known BH potential Vbh- 

This latter was introduced in [TS] for the case of a static and spherically symmetric BH. In 
our approach, this case corresponds to a vanishing x% an d to the choice of the divergenceless 
vector Hf to be the gradient of an harmonic function: 

Xl = 0, H? = d t H a . (2.19) 

In the case of spherical symmetry, H a has a single pole where the black hole horizon resides 
(i.e. in the origin): 

H a (x,y,z) = h a + P a r, h a = const, r = -=L== =. (2.20) 

y x l + y z + z l 

The constants P a fit into the symplectic vector of electromagnetic charges (p A ,q\) of the BH 
itself. Within these assumptions, the tensor V%j relates to Vbh as follows: 

% = d^rVsH, V BH = ~P J MP. (2.21) 

Moreover, Maxwell Eqs. (12.151) decouple completely from (I2.16P and (I2.17p . which get exactly 
the form presented in [T5]. Therein, it was also shown that in the case in which the action (12. ip 
describes the bosonic sector of iV = 2, D = 4 Maxwell-Einstein supergravity, the BH poten- 
tial Vbh acquires a nice geometrical interpretation in terms of special Kahler geometry (see 
e.g. [19] for a review and a list of Refs.): 

Vbh = ZZ + G a ~ a D a ZD- a Z, (2.22) 
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where Z is the N = 2 central charge function and D a Z (matter charges) denotes its Kahler- 
covariant derivative: 

Z=(P,V), D a Z=(P,D a V), (2.23) 

defined in terms of the symplectic sections V a of the flat symplectic bundle of special geome- 
try [201125. 

Interestingly, the Cartesian D = 3 tensor BH potential V%j fl 2 . 1 8 [) can be given a geomet- 
rical interpretation in the spirit of Eq. f !2.22[) . Indeed, if one constructs a D = 3 Cartesian 
generalization of the central charge and of its covariant derivatives as 

Zi = (Hi,V), Zri^faDaV), (2.24) 

then it can be computed that 

Vij = Z^Zj) + G aa Z a (iZ aj ), (2.25) 

which can be regarded as the generalization of (12.221) to generic, not necessarily axisymmetric, 
and thus possibly multi-centered, stationary solutions (at least those with flat three-dimensional 
spatial slices, cfr. (12. lip ). 

For later convenience, let us also here define the non-rotating limit (xi = 0) of Z^ and Z ai 
defined in fl2^j) : 

Z % = (H h V), Z ai = (H h D a V). (2.26) 

As it will become evident from the subsequent treatment, both Cartesian D = 3 generaliza- 
tions (12.241) and (12.261) of the N = 2 central charge Z and of its covariant derivatives D a Z will 
play an important role in the construction of first order equations of motion for multi-centered 
BHs. 



3 First Order Equation of Motion 

The first order equations of motion arise from the following Ansatz: 

d . U = e u W . } d iZ a = e u Ul Xl = 2e- U e i} (3.1) 

where Wi, IT" and ti are functions of the scalar fields and, eventually, of some auxiliary ones. 
The choice of the Ansatz (13 .ip is motivated by the fact that for single- and multi- centered 
BPS, as well for single-centered non-BPS BHs, the first order equations are of the form (13. ip . 

In the following treatment, we will conveniently consider special Kahler geometry tensors 
whose indices are "flattened" as usual: 

Va = E^ a D a V, Z^. = Ea a Z ai etc., (3.2) 

by using the Vielbein of the scalar manifold (for further details, see App. |Aj. By means of (13.21) . 
the first order equations (13.11) can be recast as follows: 

diU = e u Wi, diZ a — e u U.fEa a , x* = 2e"^. (3.3) 

Due to the Einstein equations (12.161) . the functions entering the right hand sides of Eqs. (13. 3 j) 
must satisfy the algebraic constraint 

% = WiWj + Sa-Jlffil + iiij, (3.4) 



6 



expressing the relation between Vij and the D = 3 Cartesian fake superpotential Wi, and 
generalizing the relation [10] 

V BH = W 2 + AG a7i d a WdaW } (3.5) 

to which it reduces in the limit@ (12^511 - (IZ!ZT|I . 

For later convenience, let us define the complex D = 3 Cartesian vector 

Wi = Wi + Hi, (3.6) 

in terms of which the algebraic constraint f)3.4p can be recast afl (cfr. 

V$i = W(Wj) + $aaIlffi]y (3.7) 

We anticipate that the complex vector Wj will play an important role in integrating Maxwell 
Eqs. (j235jl . 

From now on, in order to avoid overloading the formulae, we omit underlining the flat indices. 
3.1 Construction of Flows 

We are now going to determine the first order flow-governing functions Wi, 1% (or Wj) and Ilf. 

A crucial step in constructing such flow-governing functions is to use their expression in the 
BPS case, and then perform a suitable flipping of some (linear combinations of) electromagnetic 
charges. 

3.1.1 BPS 

As it is known JTOj, [T], in the BPS class [8], at spatial infinity (i.e., where one may think of a 
restoration of the spherical symmetry), the first order superpotential W is just a combination 
of the ADM mass M and the NUT charge N, while the functions IP are the scalar charges^) 

W = Re(M + iN), IT 1 = vr a . (3.8) 

In turn, at spatial infinity it holds 

M + iN = Z, 7r a = Za, (3.9) 

and thus (13 .Sp can be rewritten as: 

W = ReZ, IF = Za. (3.10) 

In order to proceed further, we will now specialize our treatment to the N = 2, D = 4 stu 
model [16] , which can also be regarded as a common sector of D = 4 supergravity theories with 
rank-3 symmetric scalar manifolds. 

Within this model, in order to restore manifest duality covariance in f)3.10p . one has to 
consider the proper action of the compact symmetry [[/(l)] 4 = x U(l), where H4 is the 

3 Within the same limit, the first order constraint (|2.12|) reduces to the Hamiltonian constraint given by 
Eq. (11) (withc = 0) of [TS]. 

For simplicity's sake, starting from below (|3.7p we will refrain from underlining the scalar flat indices. In 
presence of curved indices, we hope the distinction will be clear from the context. 

5 Hcrc, we omit the spatial index i because in the case under consideration the spherical symmetry reduces 
the number of independent components of a Cartesian vector to one. 
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stabilizer of the completely factorized D = 4 scalar manifold [SL(2, M.)/U(l)} 3 , or equivalently 
the maximal compact subgroup (mcs) of the D = 4 generalized electric-magnetic (tZ-)duality 
group Gi = [SL(2, M)] 3 , while the commuting U (1) is the mcs of the SL(2, M.) Ehlers symmetry 
determined by the reduction to tree dimensions. As discussed in [T], under the resulting [[/(l)] 4 
symmetry, the central charge and its derivatives transform as follow^] 

Z^e 2( « >Z, Z a ^e ^« Z a , (3.11) 

while the ADM mass M, NUT charge N and scalar charges 7r a transform as 



M + iiV^e 2 » (M + iJV), 7r a ->■ e V ^• / 7r a . (3.12) 

Therefore, in order to restore manifest covariance in Eq. (13.101) . one applies transformations (I3.11j! 
on the left and right hand sides of Eq. (13.91) : 



e 2 V V J(M + iN) = e 2 V « Vz, e V »^a > / 7r * = e a V (3.13) 

by some trivial algebra, and taking into account Eq. f 13 . 8 1) . one achieves the following result, 
depending only on the phase ao: 

W = Re (e- iao Z) , n a = e ia °Za- (3.14) 

Although Eqs. ( I3.14|l were obtained in [1] at spatial infinity only, we will postulate that they 
are valid not only at the spatial infinity, but all along the corresponding whole scalar flow (with 
generally broken spherical symmetry), and we will study the consequences of this approach 
in the next Sections. Our approach then justifies the following "Cartesian generalization" 
of fl3H]): 

Wi = Re (e~ iao Z^ , IT/ = e ia °t ai , (3.15) 

with Zi and Z ai defined in ( I2.24p . Consequently, the function and hence Xii is easily deduced 
from the Einstein constraint (j3.7p : 

U = Im (e- iao Zi) X i = 2e- u Im (e^zi) . (3.16) 



Therefore, one realizes that the complex vector M + iN (composed by the ADM mass M and 
the NUT charge N) can be "prolonged" all along the scalar flow by the Cartesian D = 3 
vector Wj (13. 6p . which in this case reads 

W i = W i + i£ i = e- ia °Z i . (3.17) 

The phase a becomes a dynamical field and the consistency of the first order equations require 
that it satisfy the following equation 

diOo = e u lm(e- iao Z^ - lm(d a Kd iZ a ). 



6 The phases ao and a a (a = 1,2,3) are not directly related to the Ehlers U(l) and to the three {/(l)'s in 
i?4, but rather they are a linear combinations of the corresponding phases thereof; for further detail, see [T]. 
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It is here worth stressing that Eqs. (13. 16[) resemble the ones obtained in [8], but there is 
an important difference: they are indeed expressed in terms of Zj, instead of Zi. The very 
definition (I2.24p of the D = 3 Cartesian N = 2 central charge Zi involves terms proportional 
to the vector x%i thus Eq. (I3.16P can actually be solved for Xi itself, achieving the following 
expression: 

2e- u lm[e- iao Zi\ 
Xi ~ l-2e- u lm [ e - ia °(b, V)}' ^ 3 ' 18) 

Then, the integration of Maxwell Eqs. (I2.15P exactly implies the result of [8]. In order to achieve 
this, one is hinted by the fact that expression (13.151) can be recast in the following form: 

w; = Re<# i ,e- ia no, 

and one can thus check that 

b a = 2e u Re {e~ iao V a ) (3.19) 
satisfies the Maxwell equations (12.151) . thus yielding the following flow-defining functions: 

Wi = Re (e-^Zi) , 1% = e ia ° Z- ai} X i = -2e~ u lm (e^ ^) , (3.20) 

which govern the first-order formulation of scalar flows of the BPS class of multi-centered 
BHs 0. 



3.1.2 Almost BPS 

In order to construct a first order formalism for non-BPS flows, we will exploit suitable charge 
flippings at spatial infinity, and then consistently extend them all along the flow. The relevant 
flippings of charges have been derived in PQ by exploiting an analysis of the relevant nilpotent 
orbits of 5*0(4,4), which is the D = 3 duality group of the stu model. 

We start and consider the central charge Z and its flat derivatives Z a at spatial infinity. 
Without loss of generality (by = U(l) 3 duality), the asymptotical values of the scalar 
fields z a can be set to 

z a = -i. (3.21) 

Therefore, at spatial infinity the central charge along with its flat Kahler-covariant derivatives 
respectively read 



2^ 



1 



2^2 



q - ip 



+ Y j {p h + iq b )-p a + iq a \. (3.22) 



b ^ a 



Instead of dealing with the charges {p°,p a ,qo,Qa), let us define a basis (D°, D a , D , D a ) whose 
interpretation in terms of D-brane charges is given below pQ: 



D6 : D° = 
D2 : D n ee 



m-.D a = P a , 

DO : D EE q 0] 



(3.23) 
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this is a more convenient basis on which one can act with the charge flipping. In terms of 
D-charges (13. 23p . (13.221) can be rewritten as 



Z 

z, 



1 



2y/2 
1 



a 2V2 



D -iD° + J2(D a + iD a ) 



D -iD°-J2 ( Db + iDb) + D a + iD a 



or equivalently, inverting in terms of the Z)-charges: 



D = —= Re 
1 

72 



1 



Z + ^ZJ, D a = —lm Z + Z a -J2z i 



D° = --Llm\z + J2Za], D a = -^=Re\z + Z a -J2Zb 

b^a 



V2 

1 

72 



b^a 



(3.24) 



(3.25) 



Now, consistent with the analysis of [T], let us flip the sign of the brane charge D (I3.23|) 
in (I3.24p . obtaining 



Z = Z\ 



1 



7 = 7 



D -iD° + J2iD a + iD a ) 



D -iD°-J2 ( Db + iD b) + D a + iD a 

b^a 



by plugging the expressions (13.251) into (I3.26p . one achieves the following result: 



Z =4 



3Z - Z - 2 ReJ2 z a\ , Z a = -- \ReZ + Re J] Z b - 2Z a 



(3.26) 



(3.27) 



which can then be plugged into the right hand side of Eqs. (I3.9p . By restoring manifest (U(l)) 4 
covariance in the way described above, and extending these relations along the whole flow and 
not only at spatial infinity, one obtains the following expressions, explicitly depending on all 
four phases «o and a a 's: 



Wi = ±Re 



-tao 



3-e 



§ (a +a a - J2 a b) 



Zi-ll + e 



Ee ia «Z ai 



U a = l^z ai - \e b * a Re ( e Z, 

- b * a Im 



Ee c ^ Z bi 

. b 



(3.28) 



This case corresponds to the so-called almost BPS class [3J H], whose first order formulation 
will be studied in Sec. [5j 
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3.1.3 Composite Non-BPS 

Another possibility is to flip not only D , but also D a brane charges (I3.23P as well, namely: 



D -> -D , D a -> -D a . 



(3.29) 



By applying the very same procedure described above for the sign flip (I3.29P , one finally gets 
the following expressions, explicitly depending on only three phases a a 's: 



Wi 



n? 



b a 



(3.30) 



This case corresponds to the so-called composite non-BPS class [T], whose first order formu- 
lation will be studied in Sec. HI 



4 Composite Non-BPS Class 

We will start analyzing the composite non-BPS [TJ class ( I3.30p . its structure being relatively 
simpler than the one of the almost-BPS class. 

In the previous Section, we constructed the functions Wi and IT" f l3.30j> . governing the first 
order flow Eqs. ( 13. ip . From Einstein Eqs. ( I2.16p . one can deduce the following expression for 
the remaining flow-governing function, namely the D = 3 Cartesian rotation vector l^. 



— Im 
2 



e a Zi 



e iaa Z ai 



(4.1) 



As it can be realized by looking at previous equations, the various a-phases play an important 
role in the first order formalism for the multi-centered scalar flows. Nevertheless, their own 
dynamics has not been concerned so far. In the case This can be uniquely determined if one 
imposes that after differentiation of the first order equations (13. 3p . (13.301) and (I4.ip . the second 
order ones (12.171) are obtained. This requirement yields the following three partial differential 
Eqs. for the three phases a a 's (a = 1, 2, 3) in the composite non-BPS multi-centered class^: 



diOC a 



- e u Im 
2 



e b 



2e" 



bi 



b^a 



-ia b 



- e iaa Z n 



b^a 



(4.2) 



Thus, one can conclude that Eqs. (13. 3p . ( I3.30p . ( 14. ip . and (14. 2 p - together with Maxwell 
Eqs. ( I2.15P - constitute the first order formalism for composite non-BPS multi- centered BE 
solutions in iV = 2, D = 4 stu model. 



7 It is worth noticing that both sets (|4.2|) and (|5.2p of partial differential Eqs. for the phases (respectively in 
the composite non-BPS e almost BPS classes) exhibit a manifest stu triality symmetry [16] . 
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Let us now integrate further Maxwell Eqs. (12.151) by exploiting the hint described in Sec. 
I3.1.1l for the BPS flow. Namely, one can easily notice that for the composite non-BPS class the 
complex D = 3 Cartesian vector Wj ( 13. 6 p can be rewritten as 

Wi=(Hi,T), (4.3) 

where the complex symplectic vector 

was defined. After a little algebra, one also finds that electromagnetic potential b (12. 6p enjoys 
an expression similar to the one holding for the BPS class (I3.19p . namely: 

b a = 2 e^Re [T a (l - iBe 2U )] , (4.5) 

for some function B, which can be considered as the D = 5 rotational contribution to the 
electromagnetic potential. 

As it occurred in the BPS multi-centered class treated in Sec. 13.1.11 also in the composite 
non-BPS class Eq. (14.11) contains x% i n both sides. Remarkably, in this class the dependence on 
Xi drops out when replacing (14. 5 p into (14. ip : as a consequence, Eq. (14. ip turns into the following 
algebraic constraint: 

Im (e^ a Z l - £> ia °Z ai j = 0, (4.6) 
and the D = 3 Cartesian rotation vector Xi can be computed to read 

Xi = d i B + Be u ReLe^ aa Z i -^e ia °Z a ?\ - 2e~ u lm (e ^ z)j , (4.7) 

where the function B occurred for the first time in (14. 5p . and it should be such that Xi is 
divergenceless (c/r. (I2.14p ). 

The counting of the two "flat" directions along the (non-BPS, "large") flows in the stu 
model [121 123] is retrieved by considering that in this class the tree partial differential Eqs. ( 14. 2 p 
for the tree phases a a 's are supplemented by the algebraic constraint (14.61) (see also the treat- 
ment of next Subsection). 

4.1 Single-Centered Solutions 

Consistent with pQ (see also e.g. the two-centered analysis in [22]), we will now show that the 
well known case of a single-centered non-rotating non-BPS BH can be obtained by performing 
a suitable limit in the multi-centered class under consideration. 

We consider the single-centered limit (I2.19p - (l2.20p . in which the divergenceless vector Hf 
introduced in (12.151) is realized as a derivative of an harmonic function H a defined in (I2.20p : 

Hf = diH a , dAH a = 0. (4.8) 

Before proceeding further, it is worth clarifying here the meaning of the constraint (14.61) . With 
such a purpose, let us introduce the three real values 

Im (p i/3a z a ) 1 
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and construct an stw-like symplectic section 

V x a = (1, A 1 , A 2 , A 3 , -A 1 A 2 A 3 , A 2 A 3 , A^ 3 , A X A 2 ) T . (4.10) 

This allows the following rather elegant rewriting of the constraint ()4.6|) : 

(P,V x ) = 0. (4.11) 

Note that the three A a 's satisfy the single constraint (14. 11 p . and therefore only two of them are 
independent; this reflects the existence of two "flat" directions along the non- rotating single- 
centered stu non-BPS flow [12, 23J. A noteworthy feature of A a 's is that they are in fact integrals 
of motion in D = 3 flat space; indeed, from the above Eqs. one can check that 

di\ a = 0. (4.12) 

By inverting the definitions (14. 9p . one can thus express the three phases a a 's in terms of the 
scalar fields z a and of the integrals of motion A a 's themselves: 

e b * a = . (4.13) 

A a - z a 

In the case of spherical symmetry, each of the flow-defining D = 3 Cartesian vectors Wi and 
can actually be reduced to only one independent function: 

Wi = diTW, II? = Sir n B , (4.14) 

where r has been defined in (I2.20p : therefore, the spatial vector index i can formally be ne- 
glected. 

By plugging the phases (14. 131) into the Eqs. (I3.30p . one retrieves the known expression for 
the non-BPS fake superpotential W n sps for a non-rotating single-centered non-BPS BH in stu 
model (with non- vanishing Z at the horizon) [12J; in order to achieve a complete agreement in 
notations, one should parametrize the integrals of motion as 

A" ee e "^ + / \ (4.15) 

where the a a 's are some constant (not to be confused with the phases a a 's) satisfying a a — 0. 

a 

It should be here remarked that in the non-rotating (xi = 0) single-centered non-BPS BH 
solution supported by the Kaluza-Klein (p°, g )-configuration |2U [12] 



e~ 4U = HqH^Hs-B 2 , z a = * ; (4.16) 

B = const, At = HlI t Hz ±H H a} (4.17) 

the function Xi, defined by eq. (14. 7p . vanishes. On the other hand, if the function B is not a 
constant, then it turns out to read 

7 T cos# 
B = b + J^-, 

where J is an angular momentum; this corresponds to the Rasheed-Larsen rotating 7^ 0) 
single-centered non-BPS BH |25| 126]. 
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4.2 Multi-Centered Solutions 



4.2.1 Electric Configuration 

In PQ, it was found that for the composite non-BPS class there exists a particular solution 

e~ AU = 4 V Q H x H 2 Hz - B 2 , z a = , (4.18) 

defined in terms of functions Vo,H a and K a satisfying the following equations: 

= Hf- K a Ht + \ \e ahc \ K a K b H ci , 
d t K a = H * \ V H a H? - H b H\ + \e hcd \ K b H c H di - 2H a V K b H ci ] , (4.19) 

b 56 a b ^ a ^ c 

b ^ c 

Xi = d i B-\e ohc \d i K a H b H c , 

in which Hf, Hf and H ai = diH a are components of the divergenceless symplectic vector Hf 
introduced in (12.151) . Without entering in the details e.g. of the two-centered solution, it 
should be remarked here that in the composite non-BPS class the fact that Hf and H? are not 
harmonic implies the centers to exhibit mutually non-local electric-magnetic fluxes. The same 
holds for the magnetic solution considered in Sec. 14.2.21 

One can check that this solution satisfies the composite non-BPS first order equations derived 
above if the phases a a 's are fixed as follows: 

i^a 6 = 7r-argz a . (4.20) 

The condition (14.201) is a particular solution to the partial differential equations (14.21) . and it 
actually fixes the two non-BPS "flat" directions, whose presence cannot thus be recognized 
in the explicit form (I4.18I) - (I4.19I) of this solution. Moreover, it should be noted that (I4.20j) 
yields the vanishing of the remaining components (not entering (14.191) ) of the divergenceless 
symplectic vector Hf: 

H 0l = 0. (4.21) 

Physically, this means that the graviphoton electric charge go pertaining to each BH center 
should vanish. However, in the two-centered case, it can be checked that in general the charge 
vectors pertaining to each of the two centers are mutually non-local, even if there are no solutions 
in this class in which all two-centered duality invariants (27] are independent (also cfr. the 
analysis in [22]). 

Multi-Centered Integrals of Motion Although introduced for single-centered BH solu- 
tions, the definition (14. 9p of integrals of the D = 3 motion still holds for multi-centered BHs. 
In this case, the constraint (14. (jj) can be rewritten as 

[II Ime^](^,y A ) = 0. (4.22) 
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Assuming that the phases a a are such that the product of the imaginary parts (recall def- 
inition (14. 9p ) does not vanish, one obtains that that the skew-symmetric symplectic prod- 
uct (Hi, V\) necessarily vanishes: 

(H h V x ) = H 0i + H ai \ a - 3#fA 2 A 3) + H?\ x \*\* = 0. (4.23) 

From the definition of integrals of the D = 3 motion, and especially from its inverse ( I4.13p . one 
can see that (I4.20p holds iff all integrals of motion equal to zero: A a = 0, Va = 1, 2, 3. In order 
to fulfill the constraint (14.231) . one has to require H 0i = 0, matching (14.211) . 

The product of the imaginary parts in (14.221) might vanish. This happens for example in the 
case of the so-called magnetic configuration, which analysis is performed below. 



4.2.2 Magnetic Configuration 

If all phases a a vanish, then the constraint (14.21) turns into an identity. Despite the fact that 
the definition (14.91) of integrals of motion becomes singular, a careful analysis yields that the 
constraint (I4.22p is satisfied when 

H° = 0. (4.24) 

Physically, this condition means that the graviphoton magnetic charge p° of each BH center 
should vanish. This case was considered in {2j; therein, an explicit solution in the magnetic 
configuration 

e- 4U = W H 1 H 2 H 3 -B 2 , z a = K a + — (B — ie~ 2U ) (4.25) 

was found, in which the function Vq as well as the harmonic functions K a satisfy the first order 
equations (no summation over index a): 



diV = -H 0l - K b H bl + ^eUK'K^H*, 

H " [ H a H ai - H b H u + leUH'K^H' 1 - 2H a ^ K h H\ 



d K a = 



(4.26) 



b^a b a c 

b c 



Thus, the D = 3 Cartesian rotation vector reads 

Xi = diB + \e abc \diK a H b H c . (4.27) 

This solution, which is a particular solution of the composite non-BPS first order equations, 
is written in terms of the functions H^, H a i and Hf = diH a , which are components of the 
divergenceless symplectic vector Hf introduced in (12.151) . 



5 Almost BPS Class 

For the almost BPS [3[ [1] class of flows, the first order equations are given by the general 
Eqs. ( 13. 3 p supplemented with the flow-defining functions Wi and 11° given by Eqs. (13.281) . 
Then, one has to satisfy the Einstein constraint (13.71) . yielding the following expression for the 
D = 3 Cartesian rotation vector £y. 



1 



Im 



3 + e 



i( a 0+J2 a a) 



Zi+ 1 



-«(ao+Z>a) 



E 



e taa Z ai 



(5.1) 
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As described in the Sec. HI if one requires the equivalence of the first order equations to 
their second order ancestors, the following four partial differential Eqs. for the four phases «o 
and a a 's are obtained: 



diCto = e u lm 



(e 5 - ^e iaa +e~ ia °(^e b * a - \))Z % 

a a 
-i T. a b 

'g«o a _|_ g _ia o)(i — g bjta ) — ^ e~ iab (l — e iaa ~ ia o'j 



Z a i 



(5.2) 



d t a a = -diao - \e u cos ^±£^ Im 



e 2 



(ao— Ol a ) 



i J2 a b 

(e - 1)^ 



-i X) a b 

-(e + 3)Z ai 



in which the index c must be different from both a and b. In turn, Eqs. (15. 2ft yield a useful 
relation: 



<9; 



Q'o 



a J = 2 e 17 Re 



Im 



i(a +E a.) 



-i(a +E ««) A 

e Z, 



. (5.3) 



In order to integrate the Maxwell Eqs. (12.151) . we exploit the same approach considered in 
Sec. HI the flow-defining functions Wi and xu as we ^ as the electromagnetic potential b a , can 
again be represented as in ( 14. 3 p and ( 14. 5 p in terms of symplectic vector T, which for the class 
under consideration has a slightly more complicated expression, namely: 



(5.4) 



When replacing the electromagnetic potential b expressed in terms of T (15. 4p into Eq. (15. ip . 
this latter is turned into a constrant: 



Im 



- / 3 + e ^°+?^ Zi + (l - e - l(Q0+ ? a ^ 



0. 



(5.5) 



In other words, as for the composite non-BPS class treated in Sec. HI if one tries to deduce Xi 
from ([EI]) recalling that (cfr. d2TT5|) and (|233J) ) 



Zi = Zi + Xi(b, V), Z ai = Z ai + Xi(b, V a ), 



(5.6) 



then it can be checked that the dependence of (15. ip on Xi drops out, thus yielding the con- 
straint (15.51) . 

However, differently from the composite non-BPS class, in the almost BPS class, the Maxwell 
equations turn out to be consistent only when two further constraints are satisfied, namely: 



I : a + J2 a a a = 71 or II : oto + J2 a a a = 2 arctan (e 2U B) 



(5.7) 
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thus giving rise, as in the t 3 model investigated in [TJ], to two distinct almost BPS branched, 
which we will separately analyzed in the next two Subsections. 

The counting of the two "flat" directions along the (non-BPS, "large") flows in the stu 
model [121 [23] is retrieved by considering that in this class the four partial differential Eqs. ( 15. 2 p 
for the four phases ao and a a 's are supplemented by the algebraic constraint (15 .5p and by the 
further condition I or II of (15 .7ft . 



5.1 Branch I 

In the case I of (15 .7p . the D = 3 Cartesian rotation vector reads 

X, = d t B + e u B Re [ ^ e ia «Z ai + 3e" iQ °zJ + e~ u Im [ £ e ia «Z ai - e'^Z, 



(5.8) 



and expressions for the flow-defining functions gets slightly simplified as follows: 



Wi = Re e- ia °Zi 

n a = - e ia °t ai - - e < ao+a ^ Im 
t 2 2 



-ICto i 



+ -j(<X0+OCa) Re ^J2e ia »Z i 



bi 



(5.9) 



This branch contains BPS multi-centered solutions [8] with a particular constraint. This can 
be proved by simply setting B = in (15. 8p . thus implying the functions Wi and II? to turn into 
their BPS counterparts (I3.15p . whit rotation \i acquiring the following form: 



Xi = e Im 



a 



Independently on the value of the S-field, the constraint (15. 5 p takes the form 



Im 



e a Zi 



Y,e ia "Z ia 



0. 



(5.10) 



(5.11) 



Thus, flBTTOj) and (IBTTT]) yield 

Xi = —2 e~ u Im 



2 e~ u Im 



(5.12) 



which is exactly the expression of the D = 3 rotation pertaining to the BPS multi-centered 
class [8] subjected to the additional constraint (15. lip . 

The actual meaning of the constraint (15.111) has yet to be clarified; nevertheless, we can here 
make some observations. Considering the corresponding single-centered case, and thus setting 
both B and Xi t° zero, this constraint splits in two parts, namely: 



Im 



e iaa Z ai 







and 



Im 



»E«o 
e a Zi 



0. 



(5.13) 



8 It is interesting to notice that each branch admits a single-centered limit of different type: while branch I 
admits a particular BPS BH with 13 — (c/r. Sec. 15. 1[) . branch II exhibits a non-BPS BH (as single-centered 
limit of the particular set of embedded composite non-BPS solutions discussed in Sec. 15. 2p . This is consistent 
with the fact that two-centered almost BPS solutions are characterized by a BPS and a non-BPS center p7l [22] . 
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In this case, at least a particular solution 



arg Z a , arg Z = ^ ar § z <* 



(5.14) 



of the constraint I of (15.71) can be shown to be consistent with the first order equations and, 
as a consequence, the duality invariant 23 |2HJ EE2] vanishes. In this sense, the constraint (15. lip 
can be regarded as the generalization, within the BPS multi-centered flow embedded into the 
branch I of the almost BPS class, of the constraint is = 0. 

At the moment, it is not known whether this branch contains other consistent multi-centered 
solutions different from the aforementioned constrained BPS solutions; we leave this issue for 
further investigation. 



5.2 Branch II 

In the case II of (15. 7p . the D = 3 Cartesian rotation vector reads 



Xi 



diB - Be u Re 



3e« Z^Ys^Za* 



+ e~ u \m. 



3 e a Zi 



(5.15) 



and the constraint (15. 5p acquires the following form: 
21m 



e » Zi 



-Be 2U Re 
- B 2 e iU lm 



i a. a ^ — ^ 

3e^ Zi-J^^Zai 



e^ aa Z r -Y,e^Z ai 



0. 



(5.16) 



Here we refrain from writing down the explicit expressions for the flow-defining functions Wi 
and Ilf, resulting from plugging Eq. (15. 7p (case II) into (I3.28p . because in this branch their 
form is not very illuminating and rather cumbersome. 

It is worth remarking that, among all possible solutions to Eqs. (I5.15l) - (l5.16p . there is a 
particular one that coincides with a particular solution of the composite non-BPS treated in 
Sec [5j indeed, if the .B-field is set to a constant and the rotation Xi vanishes, then it can be 
checked that the set of equations (I5.15p - (I5.16I) turns into the set of Eqs. (14.61) and (14. 7p . which 
describes the composite non-BPS class. 



An Example with Fixed "Flat" Directions An interesting example is provided by the 
choice of all phases a a 's to be equal. The consistency of Eqs. (15.21) within the case II of (15. 7p 
requires then that 



7T 



«1 = «2 = CK3 



+ arctan Be 



2U 



3tt 



a 



2 ' 2 

This latter exactly falls into the class of solutions [21 EE] 



arctan Be 2U . 



-4U 



AH°Z 1 Z 2 Z 3 - B- 



z a = K a - 



B-ie 
2H°Z n 



2U 



where H° is harmonic, and the functions Z a and K a satisfy the following equations: 
d t Z a = H ia - |e afec | \ H\ - ^K%H ^ K c , d,K a = ~ (H« - K^H ) , 



(5.17) 



(5.18) 



(5.19) 
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where the latter equation implies K a to be harmonic. In this case, the D = 3 Cartesian rotation 
vector Xi can be rewritten as 

Xr = d t B - 2H%K a Z a , (5.20) 

and the constraint (I5.5P reads 

H l0 + H.JC- + i WWfif - BJCZ. - 2K> mH ° = 0. (5.21) 

This solution can generate other solutions by applying the ^/-duality and Ehlers symmetry 
transformations, as considered in pQ. We leave to further future investigation the issue of 
existence of general multi-centered solutions of this class; for instance, in the two-centered case, 
such general solutions should exhibit all the duality-invariant polynomials of the "minimal 
degree" complete basis [27J as independent. 



6 Conclusion 

In the present investigation, we developed the most general first order formalism for multi- 
centered and/or under-rotating extremal BH solutions with flat three-dimensional base-space in 
the stu model [TBl[T2"] of N = 2, D = 4 ungauged Maxwell-Einstein supergravity. As mentioned, 
this is a universal sector of all N ^ 2-extended supergravity theories with symmetric (vector 
multiplets') scalar manifold, and which admit an uplift to D = 5 dimensions. 

Our procedure, which generalizes the one exploited for the simpler t 3 model in [H], sets the 
non-supersymmetric (composite non-BPS [1] and almost BPS [3j[T]) classes of multi-centered 
solutions on the same footing of the well known BPS class [HJ thus allowing for a unified 
framework for the study of the flow dynamics of scalar fields. 

We developed a tree-dimensional Cartesian formalism in which the effective BH poten- 
tial [18J, as well as the scalar-dependent central and matter charges and the first order "fake" 
superpotential [10J, are generalized in a not necessarily axisymmetric framework (cfr. ( I2.18p . 
f)2.24p . (I2.25P and ( I3.4p ). which thus allows to handle also (under-) rotating and multi-centered 
BH solutions. 

Then, by extending some spatially asymptotical expressions all along the flow, we derived 
systems of partial differential equations for the a-phases describing each class of multi-centered 
systems. The corresponding counting, supplemented by an algebraic constraint (and by a fur- 
ther condition in the almost-BPS class) is consistent with the existence of two "flat" directions 
along the non-rotating sigle-centered stu non-BPS flow [T2J I2"3"] . 

The consistency of the systems of partial differential equations on which the first order 
formalism is based has then been checked also by retrieving known solutions, in the single- 
centered limit (such as the Rasheed-Larsen non-BPS Kaluza-Klein BH [25|l26] in the composite 
non-BPS class) as well as in the multi-centered case (various solutions from [H [2] have been 
obtained as particular solutions). 

Consistent with its nilpotent orbit characterization [JJ as well as with the results obtained 
in [13] for the t 3 model, the almost BPS class [3j [1] exhibits the most involved structure: 
indeed, it is described by four independent phases «o and a a (a = 1, 2, 3), and in this class the 
consistency of Maxwell equations further imposes the constraint (15. 7p . whose two-fold nature 
gives rise to two sub-branches. Such a split can be traced back to the possibility to have both 
BPS and non-BPS BH centers in this class, which thus stands on a different footing with respect 
to the BPS and composite non-BPS classes, in which the BH centers are of the same type. 
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As mentioned in Sees. H] and El we leave for future investigation the issue of existence of 
completely general solutions (with all electric-magnetic charges switched on in each BH center) 
in the composite non-BPS and/or almost BPS classes. In the two-centered case, such general 
solutions would exhibit all independent duality-invariant polynomials of the complete "minimal 
degree" basis J2JJ E21 [30] . Moreover, it would be interesting to analyze the issue of the fixed 
distance among the BH centers in the non-supersymmetric systems [H [TJ [2] , also in relation to 
the existence of walls of marginal stability. 

It is here worth pointing out that we considered only multi-centered configurations with 
"large" BH centers, namely with a well-defined near-horizon geometry (and thus, scalar dy- 
namics) already at the Einstein (two-derivatives) level; at the moment, it is not clear whether 
some solutions with "small" BH centers (cfr. e.g. [SI [31]) can be obtained in this framework 
(for an analysis in the two-centered case at the level of duality invariants, see [22]). 

The possible extension of our approach to multi-centered systems with non-flat three- 
dimensional base-space is of utmost interest, as well; within this framework, one should re- 
cover solutions found e.g. in |4]. It would also be interesting to consider generalizations to 
over-rotating (single- and multi- centered) BH solutions (for recent advances on first order 
formalism, see e.g. |32j). 

Our approach could also be applied to the multi-centered interacting non-BPS solutions of 
maximal N = 8, D = 4 supergravity recently considered e.g. in [2], and the issue of truncability 
to the universal stu sector is also worth being investigated in greater detail. 
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A Notations and Useful Formulae 

We here recall some useful notations and formulae, used throughout the present investigation. 

In a generic special Kahler geometry, the symmetric real matrix M, which firstly occurred in 
our treatment in Eq. (12. 7p , is constructed from the coupling matrices /xae and z/ae as follows [HI 



29J 




and it is symplectic: 



Pi - 



with respect to the skew-symmetric symplectic metric 




which allows one to define the symplectic product of two vectors as: 



{A, B) = A a n aP B^ = A A B A - A A B A . 
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The inverse of the coupling matrix can be easily calculated to read 

/ ,,AS „AA'„ 
V -^AA'Ai A* AS + ^AA'Ai ^S'S 

The special geometry sections are normalized as 

(V, V) = -i, (V, D a V) = 0, (D a V, D- a V) = iG a - a . 
The coupling matrix M satisfy the following identity 

l -{M aP -in af} ) = -(n aa> v a ')(n^') - (n aa ,D a v a )G a ' a (n^D^'). 

The C-tensor of special Kahler geometry and the metric of the vector multiplets' scalar manifold 
enter the basic relations 

D a D b V = iC abc G c5 D s V, D a D- b V = G ab V. 

In the stu model 

F = D abc z a z b z c = stu, 
the Vielbein and its inverse (we maintain here the underlining of flat scalar indices) 



E a - = diag ( — ^— , — , J , E a a = diag (s - s,t - t,u - u) 

\s — st — tu — UJ 

satisfy the usual definition 

G a - a = EfS^Ef, 5aa = diag(l, 1, 1). 
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